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Abstract 

We propose a simple method of explicit description of families of closed geodesies on a triaxial 
ellipsoid Q that are cut out by algebraic surfaces in M"^. Such geodesies are either connected 
components of spatial elliptic curves or of rational curves. 

Our approach is based on elements of the Weierstrass-Poncare reduction theory for hyperel- 
liptic tangential covers of elliptic curves, the addition law for elliptic functions, and the Moser- 
Trubowitz isomorphism between geodesies on a quadric and stationary solutions of the KdV 
equation. For the case of 3-fold and 4-fold coverings, explicit formulas for the cutting algebraic 
surfaces are provided and some properties of the corresponding geodesies are discussed. 

1 Introduction 

One of the best known classical integrable systems is the geodesic motion on a triaxial ellipsoid 
Q <Z M?. By introducing ellipsoidal coordinates on Q, the problem was reduced to hyperel- 
liptic quadratures by Jacobi (0) and was integrated in terms of theta- functions of a genus 2 
hyperelliptic curve T by Weierstrass in 18 . 

A generic geodesic on Q is known to be quasiperiodic and oscillate between 2 symmetric 
curvature lines (caustics). 

It is of a certain interest to find conditions for a geodesic on Q to be periodic (closed) and 
to describe such geodesies explicitly. At first sight this problem has a standard solution: by 
introducing action-angle variables {/i, I2,4'It4'2\ one can define frequencies 

=f}l(/l,/2), 02 =f^2(/l,/2)- 

Then the geodesic is closed if and only if the rotation number ^2/^1 is rational. 

However, the frequencies fij are known to be linear combinations of Abelian integrals on 
the hyperelliptic curve, hence the condition on the rotation number implies a transcendental 
equation on the constants of motion and the parameters of the problem. In practice, this appears 
to be useless for an exact description of closed geodesies. 
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Such a description can be made much more explicit when the hyperelhptic curve F turns 
out to be a covering of an eUiptic curve £ and a certain holomorphic differential reduces to 
a holomorphic differential on £. Then the corresponding geodesic itself is a spatial elliptic 
curve, which covers or a rational curve and, as an algebraic subvariety in R"^ (C"^), it can be 
represented as a connected component of the intersection of the ellipsoid Q with an algebraic 
surface. In the sequel, such class of geodesies will be referred to as algebraic closed geodesies. 

Conversely, one can also show that any algebraic closed geodesic on an ellipsoid must be a 
connected component of an elliptic or a rational curve. 

Studying closed geodesies on quadrics is a classical problem. Surprisingly, we did not find 
any reference to its explicit solution in the classical or modern literature. Here we can only 
quote the paper 0, which studied the case of a 2-fold covering of an elliptic curve, when the 
solution is expressed in terms of two elliptic functions of time with different period lattices. 
Since the periods are generally incommensurable, the corresponding geodesies are not periodic 
but quasi-periodic. 

Note that in the problem of periodic orbits of the Birkhoff billiard inside an ellipsoid much 
more progress has been made (see piTlinjV 

Contents of the paper. The paper proposes a simple approach to explicit description of 
algebraic surfaces V in that cut out closed geodesies on Q. It is based on elements of the 
Weierstrass-Poncare theory of reduction of Abelian functions (see, e.g.,^1 131 U), addition law 
for elliptic functions, and the remarkable relation between geodesies on a quadric and stationary 
solutions of the KdV equation (the Moser-Trubowitz isomorphism) described in TT and recently 
revisited in j,I in connection with periodic orbits of geodesic billiards on an ellipsoid. 

Namely, for each genus 2 hyperelhptic tangential cover of an elliptic curve we construct 
a one-parameter family of plane algebraic curves (so called polhodes). Appropriate connected 
components of the polhodes have a form of Lissajou curves and describe closed geodesies in 
terms of the two ellipsoidal coordinates Ai, A2 on Q. The geodesies of one and the same family 
are tangent to the same caustic on Q and the parameters of the ellipsoid are functions of the 
moduli of the elliptic curve. 

Since equations of the polhodes depend only on the symmetric functions Ai + A2, A1A2, 
they can be rewritten in terms of Cartesian coordinates in thus giving the above mentioned 
algebraic surfaces V. 

Our family of closed geodesic contains special ones, which have mirror symmetry with respect 
to principal coordinate planes in . In particular, for the case of 3-fold and 4-fold coverings of 
an elliptic curve, the special geodesies are cut out by quadratic and, respectively, cubic surfaces 
in R'^, as illustrated in Figures [5.41 and I^TSl 

Depending on how to assign the parameters of the ellipsoid and of the caustic to the 
branch points of the hyperelhptic curve, one can obtain closed geodesies with or without self- 
intersections on Q. 

2 Linearization of the geodesic flow on the 
eUipsoid and some elUptic solutions 

We first briefly recall the integration of the geodesic motion on an n-dimensional ellipsoid 



Q 




+ 



+ 



< ai < • • • < a„ < a„+i. 
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More precisely, it is a connected component of a real part of an elliptic curve 
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Let t be the natural parameter of the geodesic and Ai, . . . , A„ be the ehipsoidal coordinates on 
Q defined by the formulas 

y2 ("i - Ai) ■ ■ ■ (flj - A„) 

Xi=ai ^ , z=l,...,n+l. 2.1 

In these coordinates and their derivatives dXk/dt the total energy 2 ("^' t^^es a Stackel form, 
and after time re-parameterization ^ 

dt = Xi--- X„ ds, (2.2) 
the evolution of A; is described by quadratures 

A^^^dAi AJ;-MA„ _ i ds for fc = 1 , 



2^-i?(Ai) 2^-i?(A„) [0 forfc = 2,...^ 

R{X) = A(A - ai) • • • (A - a„+i)(A - ci) • • • (A - c„), 



(2.3) 



where are constants of motion. 

This implies integrability of the system by the Liouville theorem. The generic invariant 
varieties of the flow are n-dimensional tori with a quasiperiodic motion. The corresponding 
geodesies are tangent to one and the same set of n — 1 confocal quadrics Qa, ■ ■ ■ , Qc„_i of the 
confocal family 

= ^^ + --- + 2±^ = llcK"+\ (2.4) 



ai — c a„+i — c 

In particular, generic geodesies on a 2-dimensional ellipsoid fill a ring bounded by caustics, the 
lines of intersection of Q with confocal hyperboloid Qd • 

The quadratures l|2.3|l involve n independent holomorphic differentials on the genus n hy- 
perelhptic curve F = {/x^ = —R{X)}, 

and give rise to the Abel-Jacobi map of the n-th symmetric product F^"' to the Jacobian variety 

of r, 



+ j wk = wfe, Pk = (Xk, ^-R[Xk)) e r, (2.6) 



Pi p, 

Po Po 

where ui, . . . , u„ are coordinates on the universal covering of Jac(r) and Pq is a fixed basepoint, 
which we choose to be the infinity point 00 on T. 

Since ui = s+const and U2^ . . . , m„ =const, the geodesic motion in the new parameterization 
is linearized on the Jacobian variety of F. 

The inversion of the map H2.6|l applied to formulas (|2.1|l leads to the following parameteri- 
zation of a generic geodesic in terms of n-dimensional theta-functions 9{wi, . . . ,Wn) associated 
to the curve F, 

X,[s)^x,-—— -, i = \,...,n + \, 2.7) 

6'[AJ(wi, . . .,Wn) 

where A = (5", (5'), rji — (j]",??^) G M^"/2R^" are certain half-integer theta-characteristics, the 
arguments wi, . . . , w„ depend linearly on ui, . . . , m„, and therefore on s, and Xi are constant 
factors depending on the moduli of F only. 

For the classical Jacobi problem (n — 2), the complete theta- functional solution was pre- 
sented in |18) . and, for arbitrary dimensions, in [TOj . whereas a complete classification of real 
geodesies on Q was made in [2]. 



^For n = 2 this re-parameterization was made by Weierstrass |18| . 
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Periodicity problem and a solution in terms of elliptic functions. As mentioned 
in Introduction, we restrict ourselves with the case when a geodesic is periodic in the complex 
parameter s, namely, double-periodic. This implies that the solution (|2.7() can be expressed in 
terms of elliptic functions of s. 

As an example, following von Braunmuhl 0, consider the geodesic problem on 2-dimensional 
quadric {n = 2) and suppose that the parameters a^, Cj in (|2.3fl are such that the curve F becomes 
birationally equivalent to the following canonical curve 

= —z{z — l){z — a){z — P){z — a(3), 

a,P being arbitrary positive constant. Then, as widely described in the literature (see, e.g., 
521E1I5)> r covers two different elliptic curves 



(l-a)(l-/3) 



with covering relations 



^.^'"l^^^if, (2.8) 



(z — a){z — f3) ' 



z =F y/a^ 



Thus, r is a 2-fold covering of £_ and £+. 

Both holomorphic differentials wi, 0^2 on T reduce to linear combinations of the holomorphic 
differentials on and £^ , namely 

dZ± z =F \fa^ 
— dz. 

W± w 

Then a linear combination of equations H2.3() for n = 2 yields 

z^ , f^' Z- , ^ 

dz + / az = —syap + const , 

w Joe w 

7^V^ f^- z+^ ^ 

dz + / dz — syap -f- const . 

w Joo w 

Inversion of these quadratures lead to solutions for Xi in terms of elliptic functions of the curves 
i±, whose arguments both depend on the time parameter s. Then, since their periods are 
generally incommensurable, the corresponding geodesies remain to be quasi-periodic. 

This observation shows that not any case of covering F to an elliptic curve results in closed 
geodesies on Q. In the next section we consider other types of coverings and obtain sufficient 
condition for a geodesic to be an elliptic curve. 



3 Hyperelliptic tangential covers and closed 
geodesies on an n-dimensional ellipsoid 

Consider a genus n compact smooth hyperelliptic surface G, whose affine part Ga C = (z, w) 
is given by equation 

R2n+iiz) being a polynomial of degree 2n -I- 1. The curve G is obtained from Ga by gluing 
the infinite point 00. Let {i7i(P), . . . , J7„(P)}, P = (z,w) E G be a basis of independent 
holomorphic differentials on G. One can also write fij{P) = (t>j{P)dT, where r is a local 
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coordinate in a neibourhood of P. Next, let A be the lattice in C" generated by 2n independent 
period vectors ■ ■ ■ , f^n)'^ ■ 

The curve G admits a canonical embedding into its Jacobian variety Jac(G) = C" [ui , . . . ,Un) / 

P^A{P) = / (I2i,...,l?„)^, (3.1) 
so that oo is mapped into the neutral point (origin) in Jac(G') and 



((/)i(oo), . . . ,0„(oo))^ 



is the tangent vector of G cJac(G) at the origin. 

Now assume that G is an iV-fold covering of an elliptic curve £, which we represent in the 
canonical Weierstrass form 

£ = W{u)f = 4p3(y) _ g^p{u) - 33 = 4(p - ei)(p - e2)(p - 63)} . (3.2) 

Here p(m) = p(u | i^,Lo') denotes the Weierstrass elliptic function with half-periods lo,lo' and 
u € C/{Zw + Zw'}. The parameters 32,33 provide moduli of the curve. 

Assume also that under the covering map t: : G ^ £ the infinite point 00 G G is mapped 
to u = 0. 

In the sequel we concentrate on hyperelliptic tangential coverings tt : G i— > when £ admits 
the following canonical embedding onto Jac(G) 

U ^ (Wi, ... ,Un) = U\J. 

That is, the images of £ and G in Jac(G) are tangent at the origin ^. The motion of tangential 
covering was introduced in [T7| in connection with elliptic solutions of the KdV equation 
(see also |31|3|[Tni). 

Namely, let 9{ui, . . . , Un) be the theta-function associated to the covering curve G and 9 be 
the theta-divisor, codimension one subvariety of Jac(r) defined by equation 6'[A](u) — 0, where 
[A] is the special theta-characteristic in the solution (|2.7|l . 

Theorem 3.1. (|13|1 For an arbitrary vector 'W G C^, the transcendental equation 

6i(Ua; + W)=0, x G C, (3.3) 

has exactly N solutions x = <?i(W), . . . ,x — (7Ar(W) (possibly, with multiplicity) . 

That is, the complex flow on Jac(G) in U-direction intersects the theta-divisor O or any its 
translate at a finite number of points. This property is exceptional: for a generic hyperelliptic 
curve G the number of such intersections is infinite. 

Note that in the local coordinates mi, . . . , u„ on Jac(G) corresponding to the standard basis 
of holomorphic differentials 

Wfe = , k = l,...,n, (3.4) 

w 

one hasU= (0, . . . , 0, 2)^. 

According to the Poincare reducibility theorem (see e.g., (21), apart from the curve £, the 
Jacobian of G contains an (n— l)-dimensional Abelian subvariety An~i. For n = 2 the subvariety 
is just another elliptic curve covered by G. 

Notice that for the case n = 2, explicit algebraic expressions of the covers and coefficients of 
hyperelliptic curves are known for iV < 8 (see |17p. 



'As follows from this definition, the 2-foId covers (12.811 are not hyperellipticallly tangential. 
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Double periodic geodesies on an ellipsoid. The algebraic geometrical property de- 
scribed by Theorem 13.11 gives a tool for a description of double-periodic geodesic flow on the 
n-dimensional quadric Q, which is linearized on the Jacobian of the hyperelliptic curve T in 
Section 2. Namely, let the genus n curves G and T are related via birational transformation of 
the form 

where (/3, 0) is a finite Weierstrass point on G and a is an arbitrary positive constant. Then the 
following theorem proved in ^ holds. 

Theorem 3.2. To any hyperelliptic tangential cover G t—^ £ such that all the Weierstrass 
points of G are real, one can associate an {n — 1) -parametric family of different closed real 
geodesies on an n-dimensional ellipsoid Q that are tangent to the same set of confocal quadrics 
Qd , • ■ • , Qc„_i • The parameters of the ellipsoid (ai) and of the quadrics (cj) are related to 
branch points of G via the transformation 



Remark. It is natural to consider a closed geodesic as a curve on Q and not as a periodic so- 
lution of the geodesic equations that depends on the initial point on the curve as on a parameter. 
That is, we disregard this parameter in the above family of closed real geodesies. 

Proof of Theorem 13.21 The transformation (|3.5|) sends the points oo and (/?, 0) on G to the 
Weierstrass points O = (0, 0) and, respectively, oo on V. Then, identifying the curves G and F, 
as well as their Jacobians, we find that the U-flow on Jac(G), which is tangent to the canonically 
embedded hyperelliptic curve F C Jac(F) at oo, is represented as the flow on Jac(F) which is 
tangent to the embedded F C Jac(F) at O, and vice versa. In the coordinates on Jac(r) 
corresponding to the basis (|2.5() . the latter flow has direction (1,0,..., 0)^ and thus coincides 
with the linearized geodesic flow on Q. 

This remarkable relation was flrst described in 11 as the Moser-Trubowitz isomorphism 
between stationary n-gap solutions of the KdV equation and generic (quasiperiodic) geodesies 
on an n-dimensional quadric. 

Next, let us flx a real constant d and the confocal quadric Qd of the family (|2.4(l such that the 
geodesies with the constants of motion ci, . . . , c„_i have a non-empty intersection with Q n Qd- 
In view of H2.1|l . when a geodesic X{s) intersects Q n Qd: one of the points Pi = {Xi,^i) on 
the curve F (without loss of generality we choose it to be P„) coincides with one of the points 
Ed± — {d, zLy/R{d)). Under the Abel-Jacobi map (|2.t)|l with Pq = oo, the condition P„ = Ed± 
defines two translates of the theta-divisor 

ed± = {^?[A](uT(?/2) = 0}c Jac(F), q = [ '^(X2i, . . . , r2„f S C". 

A geodesic is doubly-periodic if and only if it intersects Q n Qd at a finite number of complex 
points. In this case the linearized flow on Jac(F) must intersect Qd± at a flnite set of points too. 
In view of the Moser-Trubowitz isomorphism and Theorem 13. II this holds if F is a hyperelliptic 
tangential cover of an elliptic curve £. Then, under the transformation (|3.5|l with an appropriate 
P, the real Weierstrass points on F give real and positive parameters ai , Cj of the doubly-periodic 
geodesic. 

Finally, there is an {n— l)-dimensional family of elliptic curves £ in Jac(G), which is locally 
parameterized by points of their intersection with the Abelian subvariety An-i- This gives rise 
to an (n — l)-dimensional family of the doubly-periodic geodesies, q 



Remark. Since for any chosen A^-fold tangential cover G ^ £ the branch points of G are 
functions of the two moduli (72, 53, the parameters a.^, Cj are uniquely determined by them and by 
the rescaling factor a in (|3.5() . This implies that not any ellipsoid Q may have doubly-periodic 
geodesies associated with the given degree of covering as described by Theorem 13.21 One can 
show that even in the simplest case of a triaxial ellipsoid (n = 2) and = 3 or 4, for any 
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fixed positive 01,02 there exists only a finite number of possible 0,03 for which the geodesies 
are doubly-periodic^. 

Naturally, this does not exclude the existence of such geodesies for other degrees of tangential 
coverings or those obtained from a periodic flow on Jac(G') via a birational transformation 
different from l|3.5|l . or even just closed geodesies, which are not doubly-periodic. However, the 
latter, if exist, cannot be algebraic curves in view of the following property. 

Lemma 3.3. Any algebraic closed geodesic on ellipsoid Q C M" is a connected component of 
an elliptic or rational curve. 

Proof. Let a closed geodesic be a connected component of an algebraic curve C. Since the 
geodesic flow on Q is linearized on an unramified covering of Jac(r), C must be an unramified 
covering of an algebraic curve Cq C Jac(r) and and, moreover, Cq must be a one-dimensional 
Abelian subvariety. Then, if F is a regular curve and, therefore, Jac(r) is compact, Cq can be 
only elliptic. If F has singularities (when, for example, Cj — ai and the geodesic lies completely in 
hyperplane Xi — 0) and its generalized Jacobian is not compact, then Cq can be also a rational 
curve. In both cases C, as an unramified covering of Co, can be only elliptic or a reducible 
rational curve. |— | 

In the case n = 2 the algebraic closed geodesies on a triaxial ellipsoid can explicitly be 
expressed in terms of symmetric functions of the two ellipsoidal coordinates Ai, A2 on Q. As 
a result, such geodesies can be rewritten in terms of Cartesian coordinates in M^. We shall 
describe this procedure in the next section. 

4 Genus 2 hyperelliptic tangential covers, algebraic pol- 
hodes, and cutting algebraic surfaces in 

Suppose that the genus 2 hyperelliptic curve G 

5 

w^ = ^l[iz-bk) (4.1) 

fc=i 

is an A^-fold tangential covering of the elliptic curve £ in H3.2(l . Then, according to the Poincare 
reducibility theorem, G is also an A'^-fold covering of another elliptic curve 

^2 = {W^^- (4^3 - G2Z - G3) = -4(Z - Ei){Z - E2){Z - E:i)] , 

the parameters G2,G3 being functions of the moduli (?2,<?3- 

Let U € <Che uniformization parameter such that Z = p(/7), W — p'{U), and p is the 
Weierstrass function associated to the curve £2. As above, assume that the point 00 g G is 
mapped to U — 0. Then one can show that the map tt : G — > £2 is described by formulas 

Z = Z{z), W^w'^Wiz), (4.2) 

where fc is a positive odd integer number and Z{z), ^{z) are rational functions of z such that 
Z{oo) = 00. The second relation in 14.211 implies that the Weierstrass points on G are mapped 
to branch points on £2. 

Consider the canonical embedding of G to its Jacobian variety = (ui,U2)/A, 

\Joc W W J 

^Explicit algebraic conditions on 01,02,03,0 for the case of 3- and 4-fold tangential covers were 
presented in |T|. 
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The image of the embedding is the theta-divisor 8 that passes through the origin in Jac(G) and 
is tangent to vector U = (0, 1)"^. 

£2 is hfted to the Jacobian variety of G. Namely, for any point 

, P(^) } e G, one has 



The second covering tt : G 



p(.) 



dz 
w 



dZ 



(4.3) 



where k is a constant rational number depending on the degree N only. This implies that 
z-coordinates of the N points of intersection of a complex 'U2-liiie {ui =const) with G ~ G 
Jac(G) are the roots of the first equation in H4.2|l with Z = Z{Q) (see also 15 ). 
Now let Pi — {zi,wi), P2 = {zi,W2) G G and consider the full Abel-Jacobi map 



A{Pl)+A{P2) - {ui,U2f. 



(4.4) 



Assume that ui,U2 evolve according to U-flow, that is ui =const. Hence ^1,2:2 satisfy the 
equations 

zi = , Z2 = . (4.5) 



Zi - Z2 Z2 - Zi 

This imposes a relation between coordinates of Pi and P2 on G. In the generic case, the relation 
is transcendental one and the coordinates are quasiperiodic functions of time. However, if G is 
a tangential covering of an elliptic curve, then the relation becomes algebraic and can be found 
explicitly in each case of covering. Namely, let us set 



7r(Pi) 



Ul 



dZ 



7r(P2) 



dZ 

— , and 



Ul + U2. 



(4.6) 



In view of H4.3|l and the condition ui =const, the first equation in 1)4. 4|l implies [/, =const. 
Next, due to the addition theorem for elliptic functions. 



1 p{Ui) p'{Ui) 

1 p([/2) p'(C/2) 

1 p(-[/,) p'{-U,) 



0, 



or, in the integral form 



p((7,) + p(f/i) + p(f/2) = - 



p'{Ui) ~ p'{U2) 



1 2 



p(C/i) - p(C/2 



the coordinates Zi — p{Ui), Z2 = p(C^2) are subject to the constraint 

2G3 + G2 {Zi + Z2) + 4p(C/*) [Z2 - Zif - 4Z2Z1 [Zi + Z2) 



2 J4Z3 - G2Z1 - Ga 



4Z| - G2Z2 



Ga 



(4.7) 



Then, taking square of both sides, simplifying, factoring out (Zi — .^2)^, and replacing Zi,Z2 
by the expressions Z{zi), Z{z2) from (|4.2() . we arrive at generating equation 

16{Z{z2)-Z{zi)fpi 

+ [I6G3 + 8G2 (Z(zi) + Z{Z2)) - 32Z(zi)Z(z2)(2(zi) + 2(Z2))] p* 

+ I6G3 (Z(zi) + Z(z2)) + 8G2Z(zi)Z(z2) + 16Z{zi fZ{z2)^ + G^ = 0. (4.8) 



Written in terms of Zi, Z2, it defines an elliptic curve isomorphic to £2 for any p^ . 
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In terms of zi, 2:2, the generating equation gives a family of algebraic curves Tip, C (zi, 22), 
which we call polhodes. They are symmetric with respect to the diagonal zi = 22, as expected, 
and, for a generic p,, has degree 4iV^. A polhode describes an algebraic relation between z- 
coordinates of the divisor Pi ~ {zi,wi), P2 = (21,^2) on G, which holds under the U2-flow on 
Jac(G'). The parameter p* = p(f/*) plays the role of a constant phase of the flow. 

The polhodes thus can be regarded as ramified coverings of £2 and, therefore, in general, 
have genus > 1. 

Real finite asymmetric part of polhodes. Suppose that all the roots of the degree 5 
polynomial in l|4.1|l are real and set 

61 < &2 < • ■ • < ^>5- (4.9) 

Assume that the variables 21, 22 range in finite segments [bi, bj], where both wi, W2 are real and 
finite. Taking in mind applications to problems of dynamics, we also assume that these segments 
are different and 21 < 22- Then the motion of the point (21, 22) is bounded in the unique square 
domain 

S* = {62 < 21 < 63, 64 < 22 < 65}- 

Let also p* in H4.8|l be real. The part of polhode Hp, C (21, 22) that lies in S will be called the 
real asymmetric part of Hp,. At the vertices of the domain both u;i,?i'2 equal zero. Then, in 
view of equations (|4.5() . this part of the polhode is tangent to the sides of S or passes through 
some of its vertices. 

Lemma 4.1. //[/* is such that in the domain S one of the following relations holds 

p{U,)^Z{zi), or p{U,)^Z{z2), (4.10) 
then the real asymmetric part of "Hp, is empty. 

Proof. Indeed, in view of (|4.6|l . condition p* = Z{zi) implies U2 — 0. Hence, for the above 
value of 2i S [62,63], the coordinate 22 must be infinite. If the component of Tip, in S is not 
empty, then the polhode must intersect the boundary of S, which is not possible. Hence this 
component is empty. 

If the second condition in H4.10|l is satisfied, the proof goes along similar lines, q 

In view of the above lemma, we also assume that the constant parameter p* lies in a segment 
on R where neither of the conditions ()4.1U|I is satisfied, which is one of the gaps [Ea,Ef}], 
[-00, El], [£^3,00]. 

Special polhodes. If the parameter p* in H4.8|l coincides with a branch point of £2, then 
the equation of the polhode simplifies. 

In the first obvious case p* = cx) the generating equation H4.8|l reduces to Z{zi) — ^(22) = 0. 
Since Z{z) is a rational function, from here one can always factor out zi — 22. Thus, the 
connected component of the polhode in the domain S is 

K.. = (fe>^^0|. ,4.11) 

I 21-22 J 

Next, for p{U^) — E^, p'{U^) = 0, from the addition formula (|4.7|) we obtain the following 
simple equation 

p'iUi)iE^ - p(C/2)) = p'miE,, - p{Ui))). (4.12) 
^However, as seen from the structure of 14.81 . a fixed generic zi (and u\) results in 2N (complex) solutions Z2. 
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Taking squares of both sides and simplifying, we get 

4(Zi - Z2){E^ - Z2){Ea - Zi) ■ {E^{Zi + Z2) - Z1Z2 + El+ EpE^) = 0, 

(a,/3,7) = (1,2,3). 

Then we factor out the term {Zi — Z2) that leads to polhode Tioo, as well as the product 
[Ea — Z2){Ea — Zi) that leads to two lines in (zi, Z2)-plane and therefore cannot describe the 
polhode. As a result, we obtain the special generating equation 

(Z(zi) - E^)[Z(z2) - E^) ~ 2El - EpE^ - 0, (4.13) 

which defines the special polhode "He^- For a fixed generic zi, this equation has N complex 
solutions for Z2. 

Lemma 4.2. The polhodes Tioc, pass through two vertices of the domain S. 

Proof. Since 6 branch points of G are mapped to 4 branch points of £2, some different finite 
branch points of G are mapped to the same finite branch point on the elliptic curve. Thus, at 
two vertices of S, Z{zi) = Z{z2) for zi ^ Z2 and the polhode TLoc passes through these vertices. 

Next, at the vertices of S one has wi = W2 = and, in view of the second relation in H4.2|l . 
p'{Ui) = p'{U2) — 0. Hence equation H4.12|l is satisfied in all the vertices. On the other hand, 
at two of the four vertices the condition {Ea — U2){Ea — ui) = is also satisfied. Since the 
product (Ea — U2){Ea — Ml) was factored out in H4.13|l . the polhode TLec does not pass through 
the latter two vertices, hence it passes through the other two. 



Polhodes and Closed Geodesies on an Ellipsoid. Let (/3, 0) be a finite Weierstrass 
point on G. Then under the birational transformation [z, w) 1— > (A, /x) given by (|3.5f) with a = 1 
and (3 — hi (the minimal root of (|4.9|) ') the curve G passes to a genus 2 curve 

T^{^?^ -A(A - ai)(A - a2)(A - a3)(A - c)}, 

such that ai and c are positive. Thus T can be regarded as the spectral curve of the geodesic 
flow on the ellipsoid 

XI xi xl \ 

— ^ ^ = \\ , ai <a2 <a3 

ai a2 03 J 

and the corresponding variables 

Ai = ^^, A2 = ^^, (4.14) 

Zi — Oi Z2 — Oi 

are the ellipsoidal coordinates of the moving point on Q. 

In view of Theorem 13. 21 under the transformation (|3.5|l the real asymmetric part of polhode 
Tip, describes a closed geodesic on the ellipsoid Q in terms of the ellipsoidal coordinates, whereas 
the whole family of the polhodes gives a one-parametric family of such geodesies that are tangent 
to one and the same caustic on Q. 

Substituting expressions zi(Ai), 2:2 (A2) into the generating equation (|4.8|l . one obtains equa- 
tion of the geodesic in terms of symmetric functions Ei = Ai -I- A2, S2 = A1A2 of degree 2N. In 
view of relations H2.1|l for n = 2, the latter can be expressed via the Cartesian coordinates as 
follows ^ ^ 

51 = ai + 03 H (02 - ai)Xf H (02 - a3)X|, 

ai as 

52 = aia3 H (a2 - aiJX^ H (02 - ci^jX^. 

ai 0,3 

As a result, one arrives at equation of an algebraic cylinder surface Vp, of degree 4N in M'^, which 
cuts out a closed geodesic on Q. More precisely, one get a family of such surfaces parameterized 
by p*. 
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Remark. Since the equation depends on squares of Xi only, such surfaces are symmetric 
with respect to reflections Xi —Xi. Thus, the complete intersection Vp, n Q consists of a 
union of closed geodesies that are transformed to each other by these reflections. An example 
of such intersection is given in Figure 15.41 

As we shall see below, in some cases the equation admit a factorization and the cylinder Vp^ 
splits in two connected non-symmetric components. 

It should be emphasized that the method of polhodes is based on the existence of the second 
covering G i~> £2 and the addition law on £2 , so it does not admit a straightforward generalization 
to a similar description of algebraic closed geodesies on n-dimensional ellipsoids (n > 2) . Indeed, 
as mentioned in Section 3, in this case £2 is replaced by an Abelian subvariety An-i, for which 
an algebraic description is not known. 

In the sequel we consider in detail polhodes Hp, and surfaces Vp, for the 3:1 and 4:1 hyper- 
elliptic tangential covers. 

5 The 3:1 tangential cover (the Hermite case) 

In this case first indicated by Hermite see also 0] the elliptic curve £1 in (|3.2() is 

covered by the genus 2 curve 



G = = -i (4z3 - 9g2Z - 2733) {z^ - ^92)^ 



(5.1) 



The latter also covers the second elliptic curve 



£2 = {W^^- {4W^ - G2W - G3) = -4(Z - Ei)(Z - E2)(Z - E3)} , (5.2) 
G2 = ^-^{9l + 9gl), Gs^^-fg,{3gl-9l) (5.3) 

and the covering formulas l|4.2|l take the form 

Z = 1 (4z3 _ - 953) , W ^ {Az^ ~ 3g2) . (5.4) 

The roots of the polynomial in (|5.1|l are real iff 32,53 are real and 52 > 0, 2753 > 52- Then, 
assuming that ei < 62 < e^, Ei < E2 < E^, the following ordering holds 

h = -\/3g2, {&2, ^3, h} = {3ei, 3e2, 3e3}, 65 = ^3^2, 

^1 = (353 + , E2^-^(^3g3-^j3^y E3 = ^g3, (5.5) 

and 

Z{bi) = Ei, Z{h)=E2, Z{b2) = Zib3) = Z{bi) = E3. (5.6) 

Now substituting (|5.4|l into the generating equation H4.8|l and taking into account H5.3|l . we 
get following family of polhodes 

M2izi,Z2)pl+Miizi,Z2)p, + Mo{zi,Z2) = 0, (5.7) 
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where 



M2 = ^{Z2 ~ zif (952 - 4ziZ2 - ~ 4z|)2, 
Id 

Ml = — ^.92.93 - -^52 (^1 + ^2) - 2z|z3 (^zl - Z1Z2 + z|) (zi + Z2) 

+ ^.92 {4 + 4~ ZlZ^ - Z?22 + 3zfz|) + ^533! (421^2 + zf + z|) (zi + 22)^:12:2 

81 27 

- ^535^2 (^t + zl + 2ziz| + 2z3z2) + — gl (23ziZ2 + 4^^ + 4z|) (zi + Z2) 

81 9 

- — g| (2z2 - Z1Z2 + 2z|) (zi + Z2)ziZ2 + -93 (z^ + z^ + Iz^zf) , 

Mo - ^ {729gl + 4374g|ziZ2 + 52 ISSglffi + 256z^^) 

+ ^228^ ^3.92 (zi + Z2) - ^33 (zi + Z2) (21 - Z1Z2 + z|) zj'zl 

9 /2 2\44 ^^^"^22/'-i/~i 2 2\ 

- 2-92 (^1 + ^2) ^1^2 + ^^5352 (IOZ1Z2 + Zi + Z2) 

+ ^ffi (4 + 4 + 10z?zi) - ^52 (8z? - 27ziZ2 + 8z|) Z1Z2 

+ Y^ff2 (4 + 4 + ^44) 44 - ^4 (27z? - 4ziZ2 + 27z2) zfzj 

729 243 

- ^5352 (zt + 4 + 52iz3 + 5z3z2) - Y^5332 (20zi2 - 47ziZ2 + 20z2) (zi + Z2) 



81 

+ -^9392 (zt + Z2 - Z1Z2 - zfz2 + 3Z1Z2) (zi + Z2)ziZ2 


729 

- "^5352 {'^4 - ZlZ2 + 2zl) (Zi + Z2)Z1Z2. 

Reality conditions. Assume that the parameter p* is real and (zi, Z2) range in the square 
domain S, namely 

5 = {3ei < zi < 3e2, 3e3 < Z2 < v^3g2}- 

Then, from relations (|5.5|l . 1)5. 6|l . we conclude that the conditions H4.1(Jfl in Lemma f4. II cannot 
be satisfied if and only if varies in the gap (—00; Ei\. 

Limit Polhodes. When = 00, equation (|5.7() reduces to 

952 - 4ziZ2 - 4z? - 4z| = 0, (5.8) 

which can also be obtained directly from (|4.11|l . Thus Tioo defines a conic C in (zi, Z2)-plane, 
which passes through 2 vertices of the domain 5. 

Next, setting in (|4.13|) p» = i?i = — 1(333— v/Sgf) ^'^"^ taking into account expressions (|5.3|) 
results in cubic equation 

16z?z| - 12^/35! (-^1 + 4) - 3652^2^1 {zl + z|) 

+ 8l5^ziZ2 + 2732^%! (zi + Z2) + 162^3^33 - 27^^ = 0. (5.9) 
The corresponding polhode H-Ei also passes through two vertices of S. 

Examples of Polhodes in S. To illustrate the above polhodes, in H3.2|) we choose 

32 - 3, .ga = 0.2. 
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Then the roots of the polynomials | (4:Z^ — ^92^ — 27173) (z^ — 852) and 
AZ^ — G2Z — Gz are, respectively, 



(3.0, 2.693, -0.201, -2.492, -3.0), and (6.3, 0.9, -7.2). 



(5.10) 



The domain S*, where the corresponding variables Wi,?i'2 are real, is 



5* {-2.492 <zi< -0.201, 2.693 < Z2 < 3.0}. 



(5.11) 



If the parameter varies in the interval (—00; —7.2), then the real roots of the equation 



do not fit into the intervals in (|5.11(l . This means that the conditions (|4.10(l do not hold in 
domain S and the real asymmetric part of the polhode may be non-empty. 

Note that if belongs to other intervals on the real line, some of these conditions are 
necessarily satisfied, so we exclude this case from consideration. 

The graphs of equation H5.7|l in the domain (|5.11|) for two generic values of p, are given 
in Figure whereas the graphs of the special polhode TLao given by equation (|5.8|) and TLei 
given by H5.9I) are presented in Figure [^| 

As seen from Figure lOl generic polhodes in S intersect generic lines 22 =const and zi =const 
at 4 and 2 points respectively. 
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(a) p. = -8.0 (b) = - 

Figure 5.1: Generic polhodes in 5" C = 



11.0 

{Z1,Z2). 




Figure 5.2: Special polhodes in S. 



Closed geodesies related to 3:1 covering. Under the the projective transformation 
(|3.5I) with (3 = bi — — v^3g2, the branch points { — y/3g2, 3ei, 3e2, Sea, V352, oo} of the curve G 
transform to infinity, 4 positive numbers {01,02,03,0} and zero respectively. Given (72,61, the 
parameters 62, 63 of the elliptic curve are defined uniquely: 

ei V3 "^i , 

Then, assuming that oi < 02 < c < 03, we get 

1 1 

2B' 



as 
a2 



Ol 



3ei + B' 

1 _ 2 
363 + B " (663 - BY 

1 2 



362 + B (663 - B)2 



2B - 363 - VSRj , 
(25 - 363 + VSi?) , 



(5.12) 



where B — \/3g2- As a result, the four parameters oi < 02 < c < 03 are uniquely defined by 52 
(or B) and ei. 
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Now we apply the transformation H4.14|l with bi = —\/3g2 to the polhode H5.7I) . This yields 
an equation of a closed geodesic on Q written in terms of the symmetric functions Si = Ai + A2, 
E2 = A1A2 of the ellipsoidal coordinates. (In fact, one obtains a family of such geodesies 
parameterized by p*.) Then, making the substitution H4.15|l one obtains the equation of the 
cylinder cutting surface Vp, in terms of squares of the Cartesian coordinates Xi^X^. For a 
generic parameter this equation has degree 12, it is quite tedious and we do not give it here. 
However, the structure of a generic polhode in S* C (zi, Z2) and the correspondence between the 
sets {3ei, 3e2, 3e3, Y^352} and {ai,a2,c, 03} is already sufficient to give a complete qualitative 
description of the geodesic on Q. 

Namely, let TZc be a ring on Q bounded by the two connected components of the caustic 
Q r\ Qc and p = k : I ^ Q he the quotient of the numbers of complete rotations performed 
by a closed geodesies in lateral and meridional directions on the ring respectively (the rotation 
number). 

Theorem 5.1. 1). Under the assumption ai < 02 < c < 03, the geodesic corresponding to a 
generic polhode j5.7| ) or to the special polhodes is located in the ring TZc between planes 
= ±h, h < y/a^ and has rotation number p — 2 : 1. It touches the caustic Q D Qc at 2 
points and has one self-intersection. 

2). Under the assumption ai < c < 02 < a^, the geodesic is located in the ring TZc between 
planes Xi = zLh, h < ^Joi and has rotation number 1:2. It touches the caustic Q n Qc at 
4 points and has no self-intersections. 
In both cases the geodesic is either a 2-fold covering of the real asymmetric part of a generic 
polhode Tip^ or a 4-fold covering of that of the special polhodes. 

Note that the self-intersection point of the polhode does not correspond to the self-intersection 
point of the corresponding closed geodesic. 

Sketch of Proof of Theorem \b.l\ Under the projective transformation (|4.14|) . a polhode Tip, C 
(zi,Z2) is mapped to a polhode Tip, in — (Ai,A2), which is tangent to lines Xj = Oi, 
i ~ 1,2,3 and Xj — c. In view of relations 1)2. l|l . the point of tangency of Tip, to the line 
Xj — flj corresponds to the moment when the geodesic X{s) on Q crosses the plane Xi = 0, 
and the tangency to the line Aj — c corresponds to the tangency of X{s) to the caustic QOQc- 
Estimating ordering and number of the tangencies of Tip, in the cases a2 < c and c < 02, one 
arrives at the statements of the theorem. |— | 

An Example of a Generic Closed Geodesic. For the above numerical choice of 52, 33 
one gets B = 3, ei = -0.83054, 62 = -0.067069, 63 = 0.89761 and the formulas (or the 

images of the values in H5.10|l ) yield 

ai = 0.16667, 02 = 0.1776, c = 0.3579, 03 = 1.96703. 

(This means that the corresponding ellipsoid is almost "prolate"^.) Projections of the intersec- 
tion Q n Vp onto (XijXs)- and (Xi, X2)-planes for p* = — 11 are given in Figure lOI 

One can see that this intersection actually consists of four closed geodesies obtained from 
each other by reflections {Xi, X2) i-^ {^Xi, ±X2). Each geodesies has the only self-intersection 
point at X3 = and corresponds to the polhode in Figure I^TI (h) which is passed two times. 

It is natural to conjecture that the four geodesies are real parts of one and the same spatial 
elliptic curve which are obtained from each other via translations by elements of a finite order 
subgroup of the curve. 

®In all our numeric examples some of the branch points of the curve F are rather close to each other. Apparently, 
this phenomenon is unavoidable and is due to the projective transformation 14.1411 . 
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Figure 5.3: Projection of four symmetrically related closed geodesies onto (XijX^) and {Xi,X2)- 
planes for = —11. In the latter case the projection of the common caustic is indicated. 



Special Geodesic for Tioo- In the special case p* = oo the equation of the surface Vp 
simplifies drastically and admits the following factorization 

{a{Xi - + pX^ -S)- (a {X^ + jf + (3X1 - (j) = 0, (5.13) 
V2 



where 



1 
6 



g (B-6ei)((2B + 3ei)- V3i?), 
^(B + 3ei) (i?-3x/3ei) , 

--^ (R + 3^/3el) rVb, 
9a V / 

(i? + 3V3ei)^ 

STI (V3i? - (3ei + 2B)) [B - 6ei) 



(5.14) 



and, as above, B = v^3g2, R — ySffz — Sef- 

Equation H5.13|l defines a union of two elliptic cylinders in that are transformed to each 
other by mirror symmetry with respect to the plane Xi — 0. It appears that each cylinder is 
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tangent to the ellipsoid Q at a point {X2 = Xj, = 0) and cuts out a closed geodesic with the 
only self- intersection at this point. As a result, the special closed geodesic on Q related to the 
polhode H5.8|l is defined by its intersection with just a quadratic surface defined by one of the 
two factors in H5.13|l ). 

Remark. Note that due to the self-intersection, the special geodesic in R-^ (P^^) is a rational 
algebraic curve and not an elliptic one, as the intersection of two generic quadrics. It admits 
parameterization 

Xi = (i + /ii cos(2i/), X2 = /i2 sin(2i^), X3 = /i3sini/, G R, 

hi, d being certain constants^. 

On the other hand, in the phase space {X, X) the corresponding periodic solution has no 
self-intersections and represents an elliptic curve. Indeed, in view of formulas H2.1|l . the latter 
can be regarded as a 4-fold covering of the rational special polhode 1)5. 8|l . The covering has 
simple ramifications at 8 points that are projected to two vertices of the domain S and two 
vertices of the symmetric domain S' obtained by reflection with respect to the diagonal zi — Z2- 
Then, according to the Riemann-Hurwitz formula (see, e.g., (3]), the covering has genus one. 
The projection — {{X,X)} C'^ = {X} maps two different points of the elliptic solution 
to the self-intersection point on Q. 

For the above values of the parameters ai, a2, a^, c the 3D graph of the special geodesic is 
shown in Figure I5T^ 




Figure 5.4: Two 3D Images of the Special Geodesic Corresponding to 7i, 
^Here the parameter u is not a linear function of time t or the rescaling parameter s in 1)2. 2|l . 
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Special Geodesic for TCei- Applying the transformation (|4.14|) with bi = — \/3.g2 to the 
special polhode (|5.9() and making the substitution (|4.15() we arrive at a sextic surface in M'^ given 
by equation 

- 256BV?^i - (48B3/i/2 + 432Se2/i/| + 288B^eihfl) XlX^ 
+ (2304B3ei/f - \%2B^jl - mUB'^elfl) Xf 

+ (192BVi/2 + 5184Be?/i/2 - SMB^eihf2 - 2^92B^elhf2) Xjxj 

+ {-AB^'fl - lOSelfi ~ WSBelfi - SGB^ifl) 

+ (3B'V2 - 11 664e5/2 - 275^61/2 + 1296B^elf2 - 108B^elf2) X^ 

+ (8645^61/1 - 46 656Bet/i - 36B^/i + 31 104B2e?/i - 7776B^elfi) Xf 

+ (I944e^/| - 125^/1 + GiSBeffi - lUB^eifl ~ 324:B^elfi) X| 

- (192B3/iV2 + 576B2ei/2/2) XfX^ = 0, (5.15) 

where 

/i = ^9e? - + Si?%/3^ , /2 = ^27e? - - 9Bei + BRV3 + 3ReiVi^ . 

It cuts out a pair of closed geodesic on Q that are transformed to each other by mirror symmetry 
with respect to the plane X2 — 0. Both geodesies have a 3D shape similar to that in Figure 15^ 
each of them has the only self-intersection point for {Xi = X3 = 0). 

Note, however, that in contrast to quartic equation (|5.13(l . the sextic polynomial in H5.15(l 
does not admit a factorization, hence none of the above geodesies can be represented as the 
intersection of Q with a quadratic or a cubic cylinder. 

For the above choice of moduli (72,33 and the parameters ai,c the projection of the sextic 
surface and the corresponding geodesies onto (Xi, X3)-plane are given in Figure 



0.5 



0- 



-0.5 



-1 




-0.4 -0.2 0.2 0.4 



Figure 5.5: Projection of the sextic surface 1)5. 15(1 and the geodesic for Ti.Ei onto (Xi, X3)-plane. 
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Remark. As follows from the above considerations, under the condition ai < a2 < c < 
all the real closed geodesies of the one-parametric family have one self-intersection point on the 
equator {X3 = 0} C Q, and as the parameter p* ranges from —00 to Ei, this point varies from 
the Xi-axis to X2-axis. 

The case ai < c < 02 < as will be illustrated in detail elsewhere. 



6 The case of 4:1 tangential covering 

This case was originally studied by Darboux and later appeared in paper [W\ in connection with 
new elliptic solutions of the KdV equation (see also Namely, the genus 2 curve G 

4 

w^ = -{z-6ei)Y[{z-zi), (6.1) 

1=1 

with 



Zl,2 =63 + 262 ± V28e^ -I- 766263 + 406^, 



Z3,4 = 62 + 263 ± Y 286§ -I- 766263 + AOej, 

is a 4- fold cover of the curve H3.2|l . It also covers second elliptic curve 
£2 = {W^ = -4(Z - Ei){Z - E2){Z - S3)}, such that 

£■1=3 (62 - 63)^ - 6 (62 - 61) (561 - 263)^, 

E2 = -6 (62 - 63)' + 3 (62 - 61) (561 - 263)^ (6.2) 
E^ = 3 (62 - 63)' + 3 (62 - 61) (561 - 263)2 

as described by one of the formulas 



9 - 36„z - 24 {e-l + 6^) - blepi 



2 



Z = Ei+F^{z), F^^-^ ^— (6.3) 

4 z — D60, 

W = -w — 

dz 

where (a, /3, 7) is a circular permutation of (1, 2, 3). In the sequel we assume 
a = l,/? = 2,7 = 3. 

Substituting projection formulas (|6.3|) to the generating equation H4.8|) one obtains equation 
of generic polhodes of degree 16, which is much more tedious than the family (|5.7|l for the 3:1 
cover, so we do not give it here. 



The special polhodes for p^, = 00 and = Ei. Substituting projection formulas 
(|6.3|) to the special generating equation H4.11|l . we obtain algebraic equation of degree 4, 

661 {zf + zf) - 366? {zf + zl) - z\zl + (I26l(zi + Z2) -zl- z|) Z2ZX 
+ (66263 -I- 36i)ziZ2 + (546? - 6I2616263 - 288616^ - 2886163) (Zi + Z2) 

+ 9 (176263 + %e\ -f 86^) (176263 + \2e\ + Sej + 8el) = . (6.4) 
Next, substituting H6.3|) into the special generating equation (|4.13|) with Ea = Ei and taking 
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into account relation Ei + E2 + = 0, one gets the following equation of degree 8 

81 z\ - 486ei(zf z\ + z\ zf) - (3159 ej + 6804 62 eg + 3159 63^) (z| + z| zf) 

+2916 ej z\zl + 2916 ei^i (za Zi + Zj ^1) + 1458 61^2 (^2 ^1 + A. 4) 
+729 (z^ + z^) - 8748^1 (z2 zf + z^ zi) + 729 z? z^ - 437461^? {zf + z^) 



-4374 eiCi^ iz2 z{ + z^ zi) ~ 2187 ^ (zf + z^) + 16 (^104976 62" + 656100 62^^ 63 
6725025 4 2 4520529 , 3 6725025 , 4 „^ r 

62^ 63^ + 62^ 63^ + 62^ 63" + 656100 63 63^ + 104976 636 - 2 £:2 

-5 E2E3~2 Es^^ zi Z2 + ^ei(1119744 62^ + 7138368 62^ 63 + 18528264 62* 63^ 

+25021467 62^ 63^ + 18528264 62^ 63"* + 7138368 62 63''^ + 1119744 63^ + 32 £^2^ 

0,, , R , .0 2165451489 , , 
+8OE2E3 + 32£;32)(zi + Z2) + 52225560 6^6^ + 107298594 6^6j^ + — 6^6;^ 



^3^ 



+14276736 63^ 63 - 144 62 63 £^3^ - 72 63^ £2^ - 180 63^ E2 £3 - 72 63^ E3 
+ 107298594 62^ 63^ + 52225560 62^ 63^ + 14276736 626^ + 1679616 63* - 72 62^ £2^ 
-180 E2 E3 - 72 62^ £3^ - 144 62 63 E2^ - 360 62 63 E2 £3 + 1679616 62* = 0, (6.5) 
= 8 62^ + 17 62 63 + 8 e3^ 6 = 13 + 28 62 63 + 13 63^ , 

Examples of polhodes in S. To illustrate the above polhodes, in the first elliptic curve 
H3.2|) we choose 

61 = -2, 62 = -1, 63 = 3, so that El = -1728, E2 = 1152, E3 = 576, 
and the roots of the polynomial in (|6.1() become 

(-12.0, -11.649, -3.0, 13.0, 13.649). (6.6) 
As a result, the square domain S d (zi, Z2) where the corresponding variables Wi, W2 are real is 

S* = {-11.649 < zi < -3.0, 13.0 < Z2 < 13.649}. (6.7) 

If the parameter p, ranges in the interval (—00; Ei = —1728), then the conditions H4.10() do 
not hold in 5, hence the real asymmetric part of the polhode is non-empty. 

Graphs of polhodes in the domain 1)6.7(1 for a generic value of p» € (—00; —1728) is given in 
Figure IHm whereas the special polhodes defined by 1(6.4(1 and 1(6.5(1 are shown in Figure |^21 

As seen from Figure 16.11 the generic polhodes in S intersect generic lines Z2 =const and 
zi =const at 6 and 2 points respectively. 
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(a) p* = -4900 (b) p. = -7000 

Figure 6.1: Generic polhodes in C = {zi, Z2) 





(a) p. = -00 



(b) p. = -1728 



Figure 6.2: Special polhodes. 



Special Closed Geodesic for TYoo- Assuming, as above, ei < 62 < 63, we conclude that 
bi = — 6ei is the minimal root of the hyperelliptic polynomial in H6.1(l . Setting this value into 
the transformation (|4.14|l . we get 



Zl 



6eiAi + 1 



Z2 



6eiA2 + 1 



Ai ' A2 
The latter substitution transforms equation (|6.4|) of polhode Hoo to 

(1026263 - 1176? + 486^ + 486^) XjXj -Xf-Xj- A1A2 

- ISei (Ai + A2) A2A1 + 9 (176263 - 6el + Sej + Sejf XfXl = 0, 



(6.8) 



which describes the closed geodesic in ellipsoidal coordinates on Q. Next, assuming that ai < 
02 < c < 03, we find 



1 



1 



ai 



03 



661 + 62 + 263 + v/M^l^^Te^^^^TIo^' 
1 



0-2 = 



661 + 62 + 263 - y/28el + 766263 + 406^ ' 
1 



661 + 63 + 262 - x/28el + 766263 + 406^ 



661 + 63 + 262 + V28efT76e2e3 + 40e| 
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Applying formulas 14.15|l and simplifying, one obtains equation of cylinder surface Voo of degree 
6 in coordinates (XijX^), which admits the factorization 

F_(Xi,X3)F+(Xi,X3) =0, 
F± = ±h3oXf + hoaXi + h2iXfX3 ± hi^XiXl ± h^X^ + hoiX^, (6.9) 

where hij are rather complicated expressions of 62, 63, so we do not give them here. 

Thus Voo is a union of two cubic cylinders V00-, V00+ that are obtained from each other by 
mirror symmetry with respect to the plane Xi = 0. 

Example of a closed geodesic for TCoo- Under the above assumption ai < 02 < c < 03, 

the values H6.6|l lead to numbers 

ai = 3.69877 X 10~^ 02 = 0.04, c = 0.111111, 03 = 2.84981 
and the equation of the cylinder V00+ becomes 

-0.301917X5* - 0.470556XiX| + 0.652885X1^X3 + 1.87228Xi - 0.303831^3 + 0.113051X| = 0. 

(6.10) 

Its projection onto {Xi, X3)-plane and the 3D graph of the corresponding geodesic are shown in 
Figure ^ . The cylinder is tangent to the ellipsoid Q at 2 points with X2 = 0, which implies 
that the geodesic has two self-intersection points. 

Remark. The special geodesic can be regarded as a two-fold covering of the plane algebraic 
curve J-± — {F±{Xi, X3) = 0}, which is ramified at two points of transversal intersection of J-± 
with the ellipse {Xf/ai + X^/a^ = 1}. (There is no ramification at the two points of tangency 
with the ellipse.) Using explicit expressions of H6.9|l . one can show that for any value of 62,63 
the genus of equals zero. Hence, according to the Riemann-Hurvitz formula, the special 
closed geodesic is a rational curve. 

However, similarly to the special geodesic for the 3:1 covering, in the phase space {X,X) 
the periodic solution corresponding to Tioc has no self-intersections and represents a connected 
component of real part of an elliptic curve. 

General closed geodesies. For a generic parameter p* the equation of the cutting cylinder 
Vp, has degree 16 and its projection onto (Xi, X3)-plane looks quite tangled: it includes 2 
intersecting closed geodesies obtained from each other by reflections Xi 1-^ — Xi. The structure 
of the real asymmetric part of generic and the special polhodes implies the following behavior 
under the condition ai < 02 < c < 03: all the closed geodesies of the family have two centrally 
symmetric self-intersection points and as the parameter ranges from —00 to Ei, these points 
vary from the plane {X2 = 0} to {Xi = 0}. All the geodesies have rotation number p — 3 : 1. 

Similarly to the case of the 3-fold tangential covering, one can consider the ordering 
ai < c < 02 < 03, which leads to generic and special closed geodesies on Q without self- 
intersections. 



The graph was actually produced by direct numeric integration of the corresponding geodesic equation with 
initial conditions prescribed by Hfci.lUII 



22 



-0.2 -0.1 0.1 0.2 

xj , 

Figure 6.3: Projection of the cylinder Voo+ onto (Xi, X3)-plane and Special Geodesic for TL, 
Conclusion 

In this paper we proposed a simple method of explicit constructing families of algebraic closed 
geodesies on triaxial ellipsoids, which is based on properties of tangential coverings of an elliptic 
curve and the addition theorem for elliptic functions. We applied the method to the cases of 
3- and 4-fold coverings and gave concrete examples of algebraic surfaces that cut such closed 
geodesies. The latter coincide with those obtained by direct numeric integration of the geodesic 
equation. This serves as an ultimate proof of correctness of the method. Depending on how one 
chooses the caustic parameter c in the interval (ai,a3), the closed geodesies may or may not 
have self-intersections. 

Thus, our approach can be regarded as a useful application of the Weierstrass-Poincare 
reduction theory. One must only know explicit covering formulas (|4.2|l . as well as expressions 
for z-coordinates of the finite branch points of the genus 2 curve G in terms of the moduli 52 , 53 • 
To our knowledge, until now such expressions are calculated only for N < 8. 

Since the method essentially uses the algebraic addition law on the second elliptic curve 
£2 , it does not admit a straightforward generalization to similar description of algebraic closed 
geodesies on n-dimensional ellipsoids: as mentioned in Section 3, in this case £2 is replaced by 
an Abelian subvariety An-i, for which an algebraic description is not known. 
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On the other hand, one should not exclude the existence of algebraic closed geodesies on 
ellipsoids related to other type of doubly periodic solutions of the KdV equation, e.g., elliptic 
not X, but in t-variable. This is expected to be a subject of a separate study. 

Our approach can equally be applied to describe elliptic solutions of other integrable systems 
linearized on two-dimensional hyperelliptic Jacobians or their coverings. 
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